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Abstract 

A space-filling, low blockage fractal grid is used as a novel turbulence gen- 
erator in a premixed turbulent combustion experiment. In contrast to the 
power law decay of a standard turbulence grid, the downstream turbulence 
intensity of the fractal grid increases until it reaches a peak at some distance 
from the grid before it finally decays. The effective mesh size and the so- 
lidity are the same as those of a standard square mesh grid with which it 
is compared. It is found that, for the same flow rate and stoichiometry, the 
fractal generated turbulence enhances the burning rate and causes the flame 
to further increase its area. Using a flame fractal model, an attempt is made 
to highlight differences between the flames established at the two different 
turbulent fields. 

Keywords: fractal grid, premixed combustion, turbulence, high-speed 
diagnostics 



1. Introduction 

Lean premixed combustion of gaseous fuels is currently one of the most 
important technologies to achieve low pollutant emissions at high efficiencies 
in the power generation sector, e.g. [1]. Reduction of NO^ emissions is a 
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direct outcome of the lower combustion temperatures when burning a lean 
mixture and complete combustion prevents the creation of unburned hydro- 
carbons and carbon monoxide. The thermal efficiency advantage provides, 
also, the added benefit of smaller CO2 emissions. 

However, for example in gas turbines, fiame stability and fiashback are 
problematic areas, among others [2, 3]. Flame stability is related to, e.g., 
fiame extinction as the lean limit is approached and fiashback can be a prob- 
lem at low heat release rates. In both cases the turbulent fiame speed is a 
determining parameter and, in general, higher turbulent velocities are pre- 
ferred. The turbulent fiame speed is, for any given fuel, determined by the 
turbulent fiuctuations of the fiow [4], so controlling the local turbulence level 
is highly desirable in a variety of situations. In laboratory premixed fiames 
a variety of square mesh grid or perforated plate designs are used at some 
position upstream of the fiame stabilization region to generate turbulence [5]. 
Very near the grid high turbulence levels can be achieved, which, though, die 
out fast following a power law decay, so the fiame is stabilized a few mesh 
lengths downstream of the grid. 

In the present paper, we propose the use of a rather new turbulence gen- 
erator, a so called fractal grid to generate the turbulence in a premixed fiame 
experiment. Fractal grids are not only very interesting for fundamental tur- 
bulence research [6], but are potentially, also, of great practical importance 
for industrial applications [7]. Conceptually, a fractal grid differs from a 
standard grid because it creates turbulence by exciting many length scales 
of various sizes simultaneously, rather than a single one, so the underlying 
mechanism that generates turbulent fluctuations is different for both types 
of grids. The grid we use possesses a space-fllling property that relates to the 
fractal dimension, Df, of the grid having the value 2 and is elaborated upon 
in [8], where it was also shown that the space-fllling grid achieves flow ho- 
mogeneity faster than fractal grids with other dimensions Df. In the present 
work, we use a grid that was "tailored" to generate the maximum of the tur- 
bulence intensity just at the position of the flame stabilizing wire (however, 
this maximum can be generated at any practical distance by using different 
grid designs). 

The structure of the paper is as follows: section 2, gives a description of 
the design and the properties of the space-fllling fractal grid and describes 
the burner, the flames studied and the diagnostics used. The results are 
presented in section 3 and a summary, conclusions and a future outlook are 
given in section 4. 
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2. Experiment 

2.1. The fractal grid 

The fractal grid, of the same design as in [8], consists of a main geometric 
pattern that is repeated at smaller scales and, as the scale decreases, the 
number of repeated patterns increases. In the grid used in the present mea- 
surements the main pattern is a square, whose bars have length Lq = 36.8 mm 
and width to = 2.70 mm. At each successive iteration of the main pattern 
there are 4 times as many squares as in the previous iteration and the length 
and width of the grid bars change as Lj = RiLq and tj = Rlto, respectively, 
where j = 0, 1, . . . , — 1 and is the number of iterations. The present 
fractal grid has Rl = 1/2, Rt = 0.4 and A^ = 3 and is shown in Fig. 1. 




Figure 1: An image of the fractal grid used in the measurements with relevant geometric 
quantities. 

The blockage ratio, a, of the fractal grid is defined as the ratio of the area 
covered by the grid to the area of the duct and it is a = 0.22. An effective 
mesh size is defined as 



a 



(1) 



where d is the burner's side (introduced in the following subsection) and P 
is the perimeter of the fractal grid; here M^ff = 13 mm. The effective mesh 
size formula when applied to a square mesh grid returns this grid's mesh size. 
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As mentioned in the introduction, one of the defining characteristics of 
this particular design of fractal grids is the difference in the decay of the 
turbulence intensity with downstream distance from the grid, as compared 
to a standard square mesh grid. In a standard grid, after a few mesh lengths 
downstream, the turbulent fluctuations decay following a power law [9] of 



the form (n ) 



~ z 



where (n ) is the variance of the velocity fluctuations. 



z is a normalized downstream distance and angle brackets denote averaging. 
In contrast, it has been shown [8, 10] that in the fractal grid used here the 
turbulence flrst increases up to a peak value at some distance downstream 
of the grid before starting to decay. The peak of the turbulence intensity 
was found to occur at a downstream distance -Zpeak ~ 0.52;*, where the wake 
interaction length scale is 

= ^ (2) 

where, for the present grid, -Zpcak = 226 mm. To demonstrate this hot wire 
anemometry is used to measure the turbulent fluctuations. Figure 2 
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Figure 2: The decay of the turbulent fluctuations, u' /U , with downstream distance. Open 
symbols correspond to the square grid and filled symbols correspond to the fractal grid. In 
this plot the point z = corresponds to the axial position of the grid. The normalization 
of the downstream distance is with z*, as introduced in the text. The vertical arrow 
shows the position of the stabilizing rod and the horizontal arrow shows the axial extent 
of the OH LIF measurements. The conditions in this plot correspond to a bulk velocity 
U = 4.3 ™/s, slightly lower than the bulk velocity at flame conditions. 

shows the downstream evolution of the centreline turbulence intensity, 
u'/U (where U is the local mean velocity), for the fractal grid and a square 
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mesh grid having the same blockage ratio and mesh size. The square mesh 
grid follows a standard power law decay whereas the turbulent fluctuations 
in the fractal grid increase with downstream distance before following a slow 
decay. Furthermore, for downstream distances larger than ~7 mesh sizes 
from the grids, both the Taylor and integral length scales are continuously 
larger for the fractal grid generated turbulence. Similarly to [8, 10] (also for 
grids with N=3), the Taylor length scale is practically constant with down- 
stream distance, whereas the integral length scale increases very slowly. So, 
the Taylor-based and turbulent (based on the integral length scale) Reynolds 
numbers have downstream evolutions of similar form to the evolution of the 
turbulence intensity. Finally, for positions at the centreline and off-centreline, 
the power spectra of the velocity fluctuations show a broad continuous power- 
law scaling region with, as expected, larger spectral densities for the fractal 
grid turbulence. 

The observed difference in the downstream behavior of the fractal grid 
turbulence, as compared to the square mesh grid, has been explained in [10] 
by considering the interaction between the wakes generated from the fractal 
grid bars, each bar having a different size; schematically, this is shown in Fig. 
3. 

□: 
a 

~^ - - (B) 

Figure 3: A schematic explanation of the interaction of the wakes generated from the 
fractal grid. Reproduced with permissions from [10]. 

The smaller wakes reach their peak turbulence intensities closer to the 
grid. They mix and would decay if it wasn't for the next size-iteration of 
wakes which reach turbulence intensities further down and thereby help the 
turbulence to build up to higher intensities further away from the grid rather 
than decay. This process continues for as many fractal size-iterations as are 
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present on the grid with the result of obtaining, without much pressure drop, 
a high turbulence intensity with good profile homogeneity at a distance z^, 
from the grid, [10, 11]. This distance, z^, is given by Eq. 2, as explained by 
the following argument [10]. The wake width, s, scales as s ~ \/tz, where t 
and z are the bar width and the downstream distance, respectively. Taking 
into account that the largest bar (with width - placed at the furthest dis- 
tance from the centreline) will create the largest wake, which will contribute 
more to the turbulent intensity, and assuming that this wake will reach the 
centreline at s ~ Lq, we can arrive at the distance z^ introduced earlier by 
equating these scalings to get Lq ~ i/toiT- 

It has been shown in [10] that the mean profile structure of the fractal- 
generated flow does not change with mean flow velocity and, in particular, 
that is independent of the mean velocity. Increasing the number N of frac- 
tal iterations improves the homogeneity of mean profiles at distances beyond 
z^. Increasing the overall combustor size and, therefore, the overall fractal 
grid size makes it easier to use a large number of iterations and therefore 
works favorably. As shown in [8, 10] who did fluid dynamic experiments in 
two wind tunnels of different sizes, scaling up does not change Eq. 2 even 
though it obviously changes 2;*, Lq and to. 

2.2. The burner 

The experimental burner consists of 2 square ducts, each of inner side d = 
62 mm, oriented vertically upwards. The upstream duct is L = 500 mm long 
and the mixture of fuel and air is injected through 4 inlets at the bottom of 
this duct. A mixing region of 200 mm is followed by a section that comprises a 
perforated plate with 4 mm holes and a mesh screen, which are used to break 
the large scale structures formed by the 4 incoming jets. After a settling 
region of 150 mm, a second flow conditioning section is placed, with a 50 mm 
long honeycomb with 4 mm holes followed by a mesh screen, in order to 
straighten and laminarize the flow and produce a uniform flowing stream of 
fuel and air. The turbulence intensity at the position where the grids are 
placed is around 1%. 

Turbulence generating grids are placed 100 mm downstream of the final 
mesh screen, exactly at the position where the second, downstream, duct is 
connected to the upstream one. In the present experiments the length of 
the downstream duct is 200 mm. A turbulent, premixed, V-shaped flame is 
stabilized on a 2.5 mm diameter rod, which is placed 50 mm downstream of 
the exit plane of the burner. The mesh size and the solidity of the standard 
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grid n'[m/s] u'/sl /[mm] Rt I /If «[°] st[™/s] 



square 
fractal 



0.19 
0.35 



0.95 
1.75 



6.4 
7.4 



76 8.9 
145 10.3 



8 
12 



0.64 
0.96 



Table 1: The flame parameters for both the square and the fractal grid. Both flames have 
the same equivalence ratio, (f) = 0.7, hence the same laminar flame speed, sl = 0.2 ™/s, 
and flame thickness. If ~ 0.72 mm. The flame angle is a, the turbulent flame speed is 
calculated as st = U sma and U — 4.6 ™/s is the area averaged velocity. The turbulent 
Reynolds number, Rt, is based on the standard deviation of the velocity fluctuations and 
the integral length scale, I. 

square mesh grid are the same as the effective mesh size and sohdity of the 
fractal grid. A general view of the burner, the laser system and the coordinate 
system is shown in Fig. 4. 



Figure 4: A sketch of the experimental setup and the coordinate system. 
2.3. The flames 

Two lean (</> = 0.7 - </> is the equivalence ratio), CH4 flames are mea- 
sured in order to present the differences between the two grids. The air and 
fuel flow rates are the same for both flames, with a bulk mixture velocity 
U = 4.6 ™/s and, at the position of the stabilizing rod, turbulence intensities, 
u'/U, 4.1% and 7.6%, for the square and fractal grids, respectively. In case 
the mean centerline velocity is used instead of the bulk velocity to normal- 
ize the turbulence intensity (see the Results section for an explanation) the 
corresponding turbulence intensities are 4.1% and 6%, respectively. Relevant 
parameters for both flames are shown in Table 1 and the position of each 
flame in the Borghi diagram is shown in Fig. 5. 

This plot shows that for the same heat release rate and the same down- 
stream distance, the fractal grid creates more intense turbulence and changes 
the combustion regime of the flame. The fractal grid flame is placed in the 
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Figure 5: The position of the two flames in the Borghi diagram, o, square grid flame; •, 
fractal grid flame. Kas = 5^Ka, where (5 = 0.1 is the thickness of the inner layer compared 
to the flame thickness, [14]. 

turbulent region of the diagram, whereas the square grid flame is at the bor- 
derline between the wrinkled flamelets and the corrugated flamelets regimes. 
This change in the relevant position of the two flames in the regime dia- 
gram implies that turbulence will play a more prominent role in the case 
of the fractal flame. In order to investigate this, high-speed laser-induced 
fluorescence imaging of OH has been applied in both flames. 

2.4- High speed OH fluorescence imaging 

The principle of planar laser induced fluorescence measurements can be 
found in [12]. A high-speed frequency doubled Nd:YAG laser (Edgewave In- 
noslab IS8II-DE) laser is used to pump a narrowband, frequency doubled dye 
laser (Sirah Allegro) that generates around 0.16 mJ per pulse at 5 kHz repeti- 
tion rate. Owing to the rather low pulse energies at this high repetition rate, 
the laser is tuned to excite the strong OH Ql(6) transition in the (1,0) band 
of the OH A-X system near 283 nm. The beam is formed into a light sheet 
resulting in laser irradiance far below saturation giving a linear dependence 
of the signal on laser pulse energy. The fluorescence from the (1,1) and (0,0) 
bands of OH is collected between 305 to 320 nm using a high-speed CMOS 
camera (LaVision HighSpeedStar 6) coupled to a two-stage high-speed inten- 
sifled relay optics (LaVision Highspeed IRO). A WG295 Schott glass fllter 
is mounted in front of the f=105, f/2.8 UV camera lens to eliminate any 
scattered laser light and a UGll fllter is used to suppress flame luminescence 
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and PAH fluorescence. The intensifier gate width is adjusted to 100 ns and 
the full extent of the array, 1024*1024 pixels^, is used. 

The thickness of the laser sheet at the measurement location is estimated 
to be 0.2 mm and the field of view of the camera is 30*30 mm^ giving a 
magnification of 0.03 ni™/pixci. The integral time scales of the flow are 1.50 ms 
and 1.45 ms for the square and the fractal grids, respectively, so that the 
measurement duration was longer than 600 integral time scales (alternatively, 
the large eddy turnover time for the square and the fractal grid is 34 ms and 
21ms, respectively, so the measurements were, at least, 30 turn over times 
long) . 

2.5. Image processing and experimental uncertainties 

The raw OH fluorescent images were processed, using the method de- 
veloped in [13], in order to extract the instantaneous distributions of the 
progress variable. The progress variable is a non-dimensional number having 
the value at the reactants and the value 1 at the products and is defined 
either as a normalized temperature or product mass fraction [14]. The flame 
surface density (here: the area of the flame per unit volume), the flame angle 
and the flame brush thickness are also calculated from these data. 

The uncertainty in the equivalence ratio and in the mixture velocity result 
from the uncertainty in the calibrated flow meters used, which have an accu- 
racy of 2.5% of the full scale. Therefore, the uncertainties in the equivalence 
ratio and the bulk mixture velocity are ±0.025 and ±0.08 ™/s, respectively. 
The statistical uncertainty in the measured standard deviation of the velocity 
fluctuations is ±8 ■ 10-^>^/s or ±16 • lO'^^/s, for the square and the fractal 
grids, respectively. The laminar burning velocity is calculated from a third 
order polynomial fit to the data of [15], as sl = 0.2 ^/s and the thermal flame 
thickness, Ip = 0.72mm, is estimated from the data of [16]. 

3. Results 

Other than the downstream evolution of the turbulence intensity, another 
difference between the square and the present fractal grid is the nozzle exit 
velocity profiles in Fig. 6, which shows the velocity across the burner at the 
exit plane, for both grids; the error bars correspond to one standard deviation 
of the velocity fluctuations at each position. 

The velocity profiles are only presented for the central part of the burner, 
extending 40 mm across, thus excluding the boundary layers and correspond 
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- e -Square grid, x-direction 

- a -Square grid, y-direction 
^ -^Fractal grid, x-direction 

-^Fractal grid, y-direction 

-%0 -20 -10 10 20 30 
Distance across the burner [mm] 



Figure 6: The velocity profile across the burner, at the burner exit, for both grids (the error 
bars correspond to one standard deviation of the velocity fluctuations) . The conditions in 
this plot correspond to a bulk velocity U — 4.3 ™/s, slightly lower than the bulk velocity 
at flame conditions. 



to a bulk velocity U = 4.3 ^/s, which is slightly lower than the mixture bulk 
velocity at the flame conditions. In both cases, the profiles show that the 
flow is symmetric and in the square grid the flow is, also, uniform across the 
burner. In contrast, the fractal grid imposes a velocity distribution across the 
burner with higher velocities at the centreline; it should be noted, though, 
that within the current experimental setup the pressure drop is the same for 
both grids, so integrating the full velocity profile (including the boundary 
layers), gives the same area averaged flow rate for both grids. The reason 
for this inhomogeneity in the case of the fractal grid is, probably, the limited 
number of fractal iterations. For fractal grids of the same design it has been 
shown [8, 10, 11] that the mean and turbulent velocity profiles across the 
mean flow direction become more and more homogeneous as the number of 
iterations increases and are, practically, homogeneous beyond ^^peak- 

Figure 7 shows the progress variable averaged over all the images, for 
both flames. 

The distribution of the fractal grid flame leans more to the incoming gases 
than the corresponding square grid flame, implying a larger turbulent flame 
speed. In fact, the flame angles measured using the c = 0.5 contour line (c 
being the progress variable) are 8 ° and 12 ° for the square and the fractal grid 
flames, respectively. Consequently, the turbulent burning velocities (shown in 
Table 1) are different for the two flames, given that the bulk mixture velocity 
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distance from wire [mm] distance from wire [mm] 



Figure 7: The mean progress variable distributions for the square grid (left) and the fractal 
grid (right) flames. 

is the same. This should be expected in view of the difference in turbulent 
fluctuations between the two flames, but it demonstrates the effectiveness 
of the fractal grid to create more rigorous burning at a given downstream 
location. 

With increased turbulent fluctuations, compared to the laminar flame 
speed, the flame will become more and more corrugated. The square grid 
flame sits at the borderline between the wrinkled and corrugated flamelets 
regimes, whereas the fractal grid flame is in the corrugated flamelets regime, 
nearer the thin reaction zones regime. Consequently, for the fractal grid 
flame, the motion of turbulent eddies will dominate over the movement of 
the flame front with the laminar burning velocity, these two effects being 
of comparable magnitude at the square grid flame. Quantitatively this dif- 
ference is shown in Fig. 8, which plots the probability density functions 
of the curvature for the two flames (positive curvatures denote flame front 
"excursions" towards the reactants). 

The pdf are shown in both linear and semi-logarithmic plots, in order 
to focus at small and large values of the curvature. The distributions are 
symmetric and similar to each other near the mean value. The tails of both 
distributions are more pronounced than a normal distribution (with the same 
mean and standard deviation) and become increasingly asymmetric at larger 
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Figure 8: The probability density functions of curvature for both flames. The semi- 
logarithmic plot at the bottom demonstrates the large differences at large curvatures 
between the two flames, whereas the linear plot at the top shows that the square grid 
flame has a larger occurrence of small and near zero values of curvature. 
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curvature values, more so for the square grid flame. The pdf of the square 
grid flame is higher and slightly wider near the mean value and it falls off 
quicker. This is more evident in the semi-logarithmic plot that shows the 
much fatter tails of the distribution of the fractal grid flame, especially at 
large positive values of the curvature. 

The downstream dependence of the flame brush thickness, 5^, is shown 
in Fig. 9. The flame brush thickness was computed from the flame surface 
density distribution (derived from the progress variable, but not shown here, 
using a procedure outhned in [13]) by fitting the sum of two gaussian func- 
tions to the transverse profile at each downstream distance and calculating 
the flame brush thickness as the average of the standard deviations of the two 
gaussian functions (each gaussian corresponds to one "leg" of the V-flame). 

8 
7 
6 

3 
2 
1 

8.1 0.2 0.3 0.4 0.5 

z/d 

Figure 9: The downstream evolution of the flame brush thickness, 6t- 

The flame brush thickness is an indication of the average movement of 
the flame due to corrugations induced by the turbulent flow field and, in 
the present flames, it is consistently thicker for the fractal grid flame (for 
distances further downstream of the stabilizing rod) manifesting the higher 
turbulence intensities in this flame. In fact, the Taylor and integral length 
scales are, also, larger in the fractal grid flame, so they induce a larger motion 
of flame. For both flames, in the current combustion regimes, large scale 
effects, depending on turbulence intensities and integral length scales, are 
expected to be more important. Smaller scale wrinkling of the flame would 
be more pronounced as one reaches the Ka = 1 line in the combustion regime 
diagram, where the flame could thicken further due to corrugations caused 
by the smallest eddies in the flow penetrating the preheat zone of the flame. 
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However, no attempt is currently made to identify such effects, which requires 
comparing the instantaneous thickness of the preheat zone with the smallest 
flow scales. 

The dynamics of the flame are, initially, explored by examining the time 
series of the flame length, as deduced from the instantaneous flame front. 
The power spectrum is shown in Fig. 10, for both grids. 

1200 



1000 

^ 800 

Zi 

a 

g- 600 



200' 



"0 100 200 300 400 500 

Frequency [Hz] 

Figure 10: The power spectrum of the flame length time series. 

As a reference, the shedding frequency from the bars of the fractal grid is 
calculated from the Strouhal number as St = ft/U, where St is the Strouhal 
number and / is the shedding frequency. Assuming a Strouhal number ^ 
0.12, the shedding frequency from the largest bar is /o ~ 200 Hz; in spite of 
that the various peaks are generally at lower frequencies than the shedding 
frequencies of the various bars. However, the peaks in the fractal grid case 
are generally of larger magnitude, implying that the fluctuations associated 
with the flame area are, also, of larger magnitude in the fractal grid flame. 
Another potentially relevant frequency that could be related to this spectrum 
is the transverse "meeting" frequency of the wakes from the various bars, but 
this frequency wasn't estimated. 

3.1. Discussion 

In the following we, briefly, assess the effect of fractal-grid generated tur- 
bulence on the turbulent burning velocity. In the corrugated flamelet regime, 
where large flow scales are more important than small scale turbulence, we 
can write [14] 
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St At , . 

where Aj- and Al are the flame surface area in the turbulent flame and A^ 
is the flame surface area of a laminar flame. With some generality, one can 
consider that the flame surface ratio At/Al is set by the ratio of a large 
(outer) length scale L^e/ (dependent on geometry) to a small (inner) length 
scale /. Our argument in this section in only indicative, and we therefore 
adopt a simple way to model this dependence as a power law 

where D — 2 > 0. If the flame was a fractal object then the exponent D 
would be its fractal dimension as proposed by [17], but we do not need to 
rely on such an assumption here. It is reasonable to expect the ratio of / to 
the integral length scale L to be an increasing function of the ratio sl/u'. 
For the sake of argument, we assume 



— ~ — , 

L \u' J 



(5) 



where the exponent h is positive. In [18] and [19] it is argued that h = 3, 
but our argument does not need this value. From the above two relations 

The above scaling is valid locally for a flamelet and we assume that we can 
relate the flame angle calculated from Fig. 7 with the local flamelet structure 
by integrating for the downstream dependence of u' and /. We can take for 
the square grid u' ~ u'q {z / Lref)~"'^'^ and l/Lref ~ (^/-^re/)^""^^, [9], and for 
the fractal grid u' ~ MqC"^^/^* (valid for distances larger than z*, see [10]) and 
an approximately constant integral length scale / ~ [u'q is the incoming 
velocity rms). We choose h = 3 and D = 2 + 1/3, mainly for convenience and 
to recover Damkohler's scaling, without this choice affecting the qualitative 
discussion and we obtain 

^T,normal ^ ^0 (7) 
STJractal ~ u'q {z^ / Lref)"^^'^ (8) 
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where St = J ST{z)dz/ L^ef- The outer scale L^e/ may be expected to scale 
with the integral scale L, which, in turn, has been shown to scale with Lq, 
[8, 10, 11]. Hence, from Eq. 2, z^,/Lre/ ~ -Z^o/^o- These expressions attest to 
an important qualitative difference between the effect on the flame, at the 
corrugated flamelets regime, from the square mesh and fractal grids, i.e. that 
the fractal grid generated turbulence imposes, also, a length scale dependence 
on the turbulent burning velocity, which can be manipulated by modifying 
the lengths and widths of the bars of the fractal grid. 

4. Conclusions 

The presented results demonstrate the effect and the value of the fractal 
grid as a turbulence generator in premixed turbulent combustion. It is shown 
that by using the fractal grid more intense turbulence can be generated at a 
given downstream distance compared to a normal square mesh grid and, as 
a consequence, the turbulent fluctuations and the turbulent flame speed are 
increased. In fact the flame angle and the turbulent flame speed increased 
by 50% by using the fractal grid. In the future, this will prove beneficial, e.g. 
for extending the lean stability limit for a given heat release rate. Analysis 
of the flame quantities reveals that both the curvature and the flame brush 
thickness of the fractal grid flame are more pronounced, in the sense that the 
flame presents larger corrugations and more intense burning than the normal 
grid flame. 

Future Outlook. The frequency content of the flame length time series shows 
some characteristic frequencies which have, however, not been analyzed in 
detail yet. We believe that further analysis of the time resolved data will 
allow the identification of dynamic effects on flame development. For ex- 
ample. Proper Orthogonal Decomposition (POD) analysis will help identify 
these dynamics, and simultaneous high-speed PIV will provide both spatial 
flow structure information (that is practically non-existent for the present 
grid design) and correlations between turbulence and combustion quantities. 
Extinction events, which are also captured within the present measurements 
and their effect on the flame will be investigated. Finally, the adverse effect 
of the mean flow inhomogeneity will be addressed, and can be ameliorated, 
by altering the design of the fractal grids. 
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